Coupling of orthogonal diffusion modes in two-dimensional nonhomogeneous systems 
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Collective diffusion coefficient in a two-dimensional lattice gas on a nonhomogeneous substrate 
is investigated using variational approach. Particles reside at adsorption sites with different well 
depths potentials and jump randomly between them. The site blocking is the only particle-particle 
interaction mechanism. It is shown that the value of the diffusion coefficient in one lattice direction 
depends nontrivially on the rate and the character of the particle jumps in all directions. When the 
jump rate in the direction perpendicular to that along which the diffusion is observed increases, the 
collective diffusion coefficient approaches values predicted within the mean field approximation. Re- 
sults of the Monte Carlo simulations for selected systems are very well reproduced by our analytical 
results. 
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I. INTRODUCTION 

Collective or chemical diffusion of adsorbed species de- 
scribes diffusion of the local density of many particle sys- 
tem involving individual jumps from one binding site to 
another. Diffusion is an important process that controls 
many physical phenomena as building of nanostructures, 
crystal growth or rate of chemical reactions. Analysis of 
such processes quite often assumes collective diffusion as 
if particles were independent of each other. Such an ap- 
proach, although reasonable as a first approximation, can 
give quite misleading results if a more precise answer is 
needed. From a perspective of a theorist collective diffu- 
sion is a complicated many-body problem of diffusion to 
which a variety of approaches are being applied ranging 
from analytic ones based on master, Fokker-Planck, or 
Kramers equations to numerical Monte Carlo or molec- 
ular dynamics simulations. An important background is 
provided by the works of Reed and Ehrlichi, an early 
summary by Gomer— , and recent reviews by Danani et 
al£ and by Ala-Nissila et alA. Relevant analytic results 
for some generic simple models were collected by Haus 
and Kehr— and interrelations between different statisti- 
cal descriptions of these processes have been reviewed by 
Allnatt and Lidiard^. 

Our interest is the coverage dependence of the collec- 
tive diffusion coefficient in a two-dimensional kinetic lat- 
tice gas model. Most of the activity in this field has been 
dedicated to interacting gases on homogeneous substrates 
with different geometries. One of the earliest seems to 
be a linear response theory approach by Zwerger 7 which 
allowed to derive analytic expressions for the coverage 
dependent collective diffusion coefficient, D(8), for a ID 
lattice gas with NN and NNN interactions. Kreuzer and 
his collaborators, using a version of the kinetic lattice 
gas model which he developed earlier to study thermal 
desorption kinetics^, investigated D(0,T) (with T be- 
ing temperature) in a ID and 2D lattice gas on nonho- 
mogeneous substrate with NN interactions and different 
models of microscopic kinetic o 10 ! 11 . 

In this work we discuss 2D systems in which interac- 



tions are limited to the site blocking only, but the local 
site potential energy landscape changes from site to site 
within one lattice elementary cell, i.e. the underlying 
substrate is nonhomogeneous. The collective diffusion 
in such a system is a complicated and difficult for an 
analytical treatment problem. Mean field treatment of 
the collective diffusion in a Schwoebel potential is due 
to Merikoski and Yin g 12 i 13 . Series of Monte Carlo sim- 
ulation data have been reported by Masin at. alJ^ and 
theoretical mean field analysis of these results based on 
approach balancing of nonequlibrium particle fluxes were 
presented by Chvoj at ali£ We have shown recently in 
Ref. [l6| that the variational approach to collective diffu- 
sion, proposed in a series of earlier work a 17 ' 18 i 19 ' 20 i 21 ' 22 , 
provides a very efficient and systematic method of ana- 
lyzing diffusion in nonhomogeneous ID systems. In this 
work we extend the results of Ref, 16 to two-dimensional 
systems. It appears that the diffusion is not a simple 
product of one-dimensional projections in two main di- 
rections of the lattice. Whereas the diffusion coefficient 
of a single particle over such a lattice always factorizes so 
the diffusion coefficients in both directions can be calcu- 
lated independently, in a many particle system the site 
blocking induces dynamical correlations between jumps 
in different directions. Properly selected variational pa- 
rameters allow to obtain the expression for the collec- 
tive diffusion coefficient which contains all possible jump 
rates present in the model. When the rate of jumps in 
one direction increases from zero to infinity, the formula 
describing diffusion coefficient evolves from the one char- 
acteristic for a one-dimensional system^ to the 2D mean 
field theory result, known from Refs. [l2lfl3l We show 
that our variational approach works very well in the de- 
scribed cases by comparing analytical results with the 
Monte Carlo data. 

Correlations of the diffusion in two orthogonal di- 
rections appear to be a very important factor re- 
sponsible to a large degree for a difference between 
the dynamical properties in narrow channels and bulk 
material o 23 i 24 ' 25 i 26 ' 27 i 28 . Transport of molecules though 
molecular pores is an essential , for its biological and tech- 



nological applications, collective process in which corre- 
lations in diffusional modes are responsible for the net 
rate of molecular transport. Results which are presented 
here allow to explore how such correlations build up. 

II. MODEL 

A system of N particles forming an adsorbate is dis- 
tributed over a two-dimensional nonhomogeneous sub- 
strate. We treat diffusion within the adsorbate using a ki- 
netic lattice gas model. Basic assumptions are standard: 
kinetics of the microstates of the lattice gas is due to 
the stochastic hopping of particles to neighbouring sites, 
only one particle in the gas hops at any given instant, 
an average residence time of particles at the adsorption 
sites is much longer than the transit time between the 
sites, the transition rates of these hops depend on the 
potential energy landscape experienced by the hopping 
particle. Double occupancy is forbidden, particles jump 
between neighboring sites with transition rates that de- 
pend on the initial and final states. Time evolution of 
this system is governed by the set of Markovian master 
rate equations for the probabilities P({c},t) that a mi- 
croscopic microstate {c} of a lattice gas occurs at time 
t 

j t p (M,t) (l) 

= E ^( W' W}) P (W], t) - W({c'}, {c})F({c}, t)\ . 
{d} 

{c} is understood as a set of variables specifying which 
particular sites in the lattice are occupied and which are 
not. W({c}, {c'}) is a transition probability per unit time 
(transition rate) that the microstate {c'} changes into {c} 
due to a jump of a particle from an occupied site to an 
unoccupied neighboring site. The rates W satisfy the 
detailed balance conditions: 

W({c}, {c'» P cq ({c'}) = W({c'}, {c}) F cq ({c}). (2) 

Here, P cq ({c}) is the equilibrium probability of a config- 
uration {c}. In the absence of interparticle interactions 
the rate depends only on the local potential energy land- 
scape experienced by the hopping particle. For thermally 
activated jumps it depends on the difference between the 
potential energy of the particle at the top of the potential 
energy barrier between the sites involved and that at the 
initial site. 

In order to investigate how the collective diffusion co- 
efficient in a given direction is controlled by the particle 
jumps and the geometry of the lattice in the direction 
perpendicular to it we analyze diffusion over two types of 
two-dimensional lattices, shown in Fig. 1. They consist 
of periodically repeated patterns of site potentials and 
intersite barriers. Elementary cell of the striped lattice 
(Fig la) is 2 x 1: there are two nonequivalent sites along 
x-axis, and all sites are identical along y-axis. Checkered 
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FIG. 1: Examples of potential geometry of the studied sys- 
tems a) striped lattice , b) checkered lattice. 

lattice (Fig. lb) has larger, 2x2 elementary cell with two 
pairs of identical sites. In general we have n x x n y unit 
cell and the potential at each site in the unit cell can be 
different. Potential energy at the tops of all barriers be- 
tween the sites are assumed to be the same. We assume 
that system under study is a box of L x x L y unit cells, 
i.e. it contains L x n x x L y n y sites. We employ periodic 
boundary conditions and eventually take a L x ,L y — > oo 
limit. 

We employ the variational approach to extract the dif- 
fusion coefficient from Eqs flTJ>I6 1 iLIS 1 i2 1 22 1 2!. To this end, 
microscopic states of the systems need to be properly pa- 
rameterized. Following Ref. [l?] we identify a microstate 
{c} by selecting one particle as a reference particle and 
specify positions of all remaining N — 1 particles with 
respect to it. Position X + alo of the reference parti- 
cle in systems with nonequivalent sites is specified using 
two variables 16 : (i) position X of the elementary cell in 
which it resides and (ii) position alo within the cell. For 
two dimensional systems X and Iq are vectors with two 
components, a is a distance between the adsorption sites 
referred to as a fractional lattice constant in what follows. 
Then a microstate {c} may be identified by the following 
set of 2(7V + 1) numbers 

{c} = [X; £ : mi, rh 2 , . . . , mjv-i] = [X; {m}], (3) 

where fhj = (m x ,m y )j is a pair of integers indicat- 
ing how far, in units of the fractional lattice constant 
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a, the j-th particle (j = 1, 2, . . . , N — 1) is away from 
the reference particle. The set of 2N integers, {m} = 
[^o : iTii, rfi2 1 ■ ■ ■ , toat_i], is referred to as a configuration, 
which accounts for the relative arrangement of particles 
in a given microstate {c}. The transition rate between 
two microstates depends on their configurations only i.e. 
WK({c}, {c'}) = W^m},^'}- This allows to take an advan- 
tage of the lattice periodicity by taking a lattice Fourier 
transform 

P {m] (k,t) = Y,e^P {rh} {X,t) (4) 
x 

of both sides of the rate equations (fTJ) . Pj^\(X,t) 
stands here for P ^{c} = [X;{m}],tV It is convenient 

to treat P^^{k,t) as an {m}-th component of a one- 
column array P(fc, t) with a macroscopic number of com- 
ponents - each component corresponds to an admissible 
microscopic configuration of the system. The Fourier- 
transformed rate equations can be written in a compact 
form 

^P(M)=M(fc)-P(M), (5) 

where "•" denotes multiplication following usual "rows 
times columns" multiplication rules. The matrix el- 
ements of M(fc) (referred to as the rate matrix from 
now on) are, in general, expressed in terms of the rates 
W{m},{m'}j except for {m} to {m'} transitions involving 
a jump of the reference atom across a boundary between 
neighboring elementary cells. For such jumps elements 
of M are multiplied by an additional k dependent factor 
exp(±fc x n x ) or exp(±fc a n a ). For details of the matrix 
M prope rties and the derivation of all formulas see Refs. 

diBQHHill. 

Eigenvalues of the rate matrix which are always neg- 
ative account for the temporal decay of a k— th Fourier- 
component of a density fluctuation from equilibrium. 
The one vanishing like \k\ 2 in the long wavelength limit, 

— Ad(/c), is referred to as diffusive eigenvalue and yields 
the collective diffusion coefficient. The corresponding 
eigenvector of M(fc) is referred to as the diffusive eigen- 
vector. This eigenvector will be calculated on using vari- 
ational formula 16 

Xl?(k)= +- [ -ff-* >X D (k) = -D\k\> ) (6) 

where <p is trial left eigenvector (possibly fc-dependent) 
and 4> is its right eigenvector counterpart with compo- 
nents 



component equal to a sum of phase factors associated 
with all occupied sites in the configuration {m}: 

N-l _ _ 

^ { ™ } (fc) = e i&( V*4) + ]T e^-^+V^V (8) 

i=i 

The phase contributed by the j-th particle is determined 
not only by its distance arhj from the reference particle 
(it is a sole contribution to the phase for a homogeneous 
system). It receives two additional distinct contributions 
Sf = (8x?.,Syi.) an d Ag = (A X £ . , A y f .) which play a 
role of the variational parameters allowing to minimize 
A^f' (fc). Both depend on the position aij within an ele- 
mentary cell of the site at which the j-th particle resides. 
The first one, 5? , called the geometrical phase, accounts 
for a nonhomogenity at the substrate within a unit cell 
and does not depend on the presence of other particles in 
the system. It is always possible to select one particular 
site (the same within each unit cell) which, if occupied, 
contributes the geometrical phase So = 0. The other 
phase, Af , called the correlational phase, is introduced 

to account for correlations between the j-th particle and 
all the remaining ones and, in principle, it depends on 
the state of occupation of all sites in the system. Fol- 
lowing previous worki 6 - we assume that it is sensitive to 
the occupation of sites nearest to lj only. Thus correla- 
tional phases A.- , associated with a pair of particles at 

£ X j — l,£ y j and ij = l X j,£yj will appear in all equations 
as a sum of contributions to the phase of the particle 
at the position ij due to its left neighbor and of the 
phase of the particle at the position £ x j — l,£ y j due to 
its right neighbor A^* = Af^._ l + A\^.^.. Similarly 

the phase A? appears with neighbors along direction y. 
Different phases associated with particle pairs are addi- 
tive. Summarizing: a phase related to each particle in 
the system depends on (i) the distance of the particles 
from the reference particle, (ii) an address within the el- 
ementary cell of the site which it occupies (geometrical 
phase), and (iii) state of occupation of the sites adjacent 
to it (correlational phase). 

Periodic boundary conditions in two dimensions imply 
in the wave number (fc) domain that the conditions 

^ik^an^L^ _ ^ ^ik y an y L y _ j 

must be used in the calculations before the long wave- 
length limit |fca| 2 <C 1 is applied. (L x - number of unit 
cells along x direction, n x - number of sites within the cell 
along x). 

We see from Eq. © that the diffusion coefficient D var 
is a ratio 



For a non-homogeneous substrate, we propose follow- 
ing Ref. [l6| that the trial left eigenvector has {m}-th 
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of the "expectation value" numerator 

no rep 
{m} ,{m'} 

x \4>\ m , } (k)-p {m} (k) 

to the "normalization" denominator 



{m} 



4>{rn}(k) 



(11) 



(12) 



Eqs. ([7]) and ([8]) have been used to get the final expression 
for the numerator in Eq. (fTTj) . Due to the detailed bal- 
ance condition each term in ([7]) correspons to a pair 
of configurations ({m},{m'}) for transitions from {rh'} 
to {m} and back. Each such pair should then appear 
in the sum only once [as indicated by the comment "no 
rep" above the sum in Eq. (flTj) ] in order to avoid dou- 
ble counting. In fact, it was shown in Ref. [H that the 
dependence of the diffusion denominator J\f(k) on varia- 
tional parameters di j and Ai j can be ignored in the long 
wavelength limit and that 



lim j\f<k) = 

k->0 



N 



/ gQVfcgT) 
^ dlnO 



= (N 2 ) - (N) 2 (13) 



Here /i is the chemical potential. The diffusion denomina- 
tor reduces to the square of the particle number fluctua- 
tion in the system, whereas the numerator M. (k) depends 
on the details of particle dynamics and on all variational 
parameters of the model. 



III. STRIPED LATTICE 

Let us consider first the striped lattice shown in Fig la. 
It consists of rows of sites, with alternating deep and shal- 
low potential energy minima. Transition rates of jumps 
along x-aie Wd for a jump from deep well and W s for a 
jump from shallow well, whereas all jumps along y-axis 
have the same rate V. A ratio r = W s /Wd between both 
rates along x-axis is the only parameter which deter- 
mines the equilibrium properties of the system at a given 
density of particles. With the site blocking interactions, 
preventing double site occupation, the mean equilibrium 
densities are 



for shallow and 



r exp(/3/i) 
1 + r exp(/3/i) 



exp(/3n) 
1 + exp(/3/i) 



(14) 



(15) 



for deep sites. The chemical potential /i controls the total 
particle density, understood as a fraction of adsorption 
sites that are occupied 9 = (9 s + 9d)/2. The denominator 




FIG. 2: Diffusion coefficient D x dependence on the total den- 
sity of the system 6 on striped (dashed line) and checkered 
lattices (full line). Jump rates along y- direction on striped 
lattices are related to the jump rates on checkered lattices for 
the corresponding diffusion curves like V = VVdVs = V s \/r. 
The lowest line corresponds to ID case, V = V 3 — Vd = 0. 
The topmost line for V = V B = Vd — > oo reproduces the mean 
field result from Refs[[lGJ. 



V=W ri r" 1 ' 2 



FIG. 3: Diffusion coefficient D x dependence on the total den- 
sity of the system on striped lattice for different rate r val- 
ues. The jump rate along y-direction is equal to the geomet- 
rical mean value of rates in x direction V = Wd/y/r. 



Af(0) given by Eq (fTJ]) in the limit (|fc|a) 2 -> depends 
only on the equilibrium properties of the system. With 
no interactions other than the site blocking we have 



(16) 



In order to derive the numerator A4(k) (fTTj) we use 
variational vector <p{ m y(k) Eq (J8]). There are two geo- 
metrical phases: <5f = (by choice) and Sf , which for 
the substrate with potential energies of all barriers be- 
ing the same is also equal to zeroi^. In the basic cell of 
striped lattice, there is only one pair of different sites, 
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which means that there is only one nonzero variational 
parameter Af . It is the correlational phase of a par- 
ticle being in a site with a shallow potential well, that 
has neighbor at left side, in deep well. All other occupa- s 
tional phases are equal to ±Af , depending on the order 
in which particles occupy deep and shallow wells. Thus 
components of trial eigenvector are ^ 

Q 


r=10- 6 

/ ,»•"'"" r-0.01 ""'\ \ 


v s =w, 


JV-l 3 

(t>{m 3 }(k) = e lk * a{m *i +A * )+lk « am i" , (17) 

2 






with Ay = 0, £f = 5 d = 0. After solving variational 
equation we get 1 


r=1 







A 6 = 



W s -W d 

w s + w d + v 



(18) 



Final expression for the diffusion coefficient along the 
x-direction in the striped system is given by 



D , 



2W S 



W d 



we s (i-e d ) 



w s + w d + w 

while for the y-direction it is 



Dy = a z V. 



,{i-e s ) + e d {i-e d ) 

(19) 
(20) 



Now, when we set V = in Eq. (JTHJ) , it simplifies 
reproducing the result for diffusion in one-dimensional 
(ID) system^ 



2W s W d 

w s + w d 



(21) 



In this limit diffusion does not depend on the density of 
the system, what is illustrated in the lowest line in Fig. 
2. For all nonzero values of V, diffusion along x-axis de- 
pends on the total density 9. This dependence changes 
with the jump rate along vertical direction V. Whereas 
boundary values that represent diffusion of the single par- 
ticle 9 = and of the single hole 9 = 1 stay unchanged 
irrespectively of the rate in ^/-direction, the height of the 
diffusion maximum grows up with the increasing value of 
V. In the limit V — > oo we have the dependence 



n inf 2 

L> = a 



2W s 9 s (l-9 d 



9 s (l-9 s )+9 d (l-9 d ) 



(22) 



This relation above reproduces exactly the for- 
mula for the diffusion coefficient in the mean field 
approximatio n 12 ' 13 . 

We can see that the rate of the jumps in the vertical 
direction V controls the character of the density depen- 
dence of the diffusion coefficient. With increasing V we 
observe a smooth transition from purely one-dimensional 
to the mean field behavior of the system. While the for- 
mer limit is obvious, the latter one can be understood in 
such a way, that a particle, capable of fast travel along 
y-axis, detects mean field occupation value of the neigh- 
boring site and instantaneously adjusts to it. 



FIG. 4: Diffusion coefficient D x dependence on the total den- 
sity of the system 9 on checkered lattice for different rate r 
values. Jump rates along y-direction are equal to the corre- 
sponding rates along x. 



Character of the density dependence of the diffusion 
coefficient depends strongly on the ratio r of rates from 
the deep and shallow sites along the direction under 
study. In Fig. 3 we show how the diffusion changes with 
r. The jump rates along y-axis increase with decreasing 
r like V = Wa/y/r, slower than the quicker of two rates 
along x: W s = W d /r. We see in Fig. 3 that with such 
a choice of parameters, the curves become higher and 
steeper as r decreases, approaching to the limiting be- 
havior in which diffusion coefficient has value D = 2W d 
for all densities, except at a one discrete point of 9 = 0.5 
for which D = 6W d . 



IV. 



CHECKERED LATTICE 



Let us now consider a checkered substrate lattice which 
has the same structure in both x- and y-directions. We 
assume, that like in chessboard every second site is dif- 
ferent (see Fig lb). Jumps out of the shallow sites in 
x-direction are given by W s , and in y-directions by V s , 
whereas jumps out of deep sites are W d and V d , respec- 
tively. The condition 



w 1 = v 1 
w d v d 



(23) 



has to be fulfilled. The equilibrium occupancies are again 
given by Eq. (TT4"|) for the shallow site and by Eq. (fT5|) for 
the deep site. Consequently, the diffusion denominator is 
again given by Eq. (fl6|) . As before, all geometrical phases 
S are equal to zero due to equal height of all potential 
energy barriers. There are now two nonzero, correlational 
phases. Minimizing the diffusive eigenvalue results in the 
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phase associated with a pair along x-direction 

X 

jw s ~w d )e s (i-e d ) 

(W s + W d )6 s (l - 9 d ) + 2V d [9 d (l - 9 d ) + 8,(1 



(24) 
(25) 



and the second phase associated with similar pair of par- 
ticles along y-direction 



A 



c 



(26) 



(V s - V d )6 s (l-6 d ) 



(y s + v d )o s (i - e d ) + 2w d [e d {\ - e d ) + e s (i - e.)] ' 

The resulting diffusion coefficient along a;-direction is 

(27) 



D c = 



2W s s (l-6 d )(W d + 2V d ) 



(w d + w s )e s (i - e d ) + 2v d [e d {\ - e d ) + e s {\ - e s )\ 

and for diffusion along y we must replace all V and W 
rates with W and V respectively. 

D x as a function of 9 for several values of rates V s 
and V d is plotted in Fig. 2. Maxima of the diffusion 
coefficient for the striped and checkered lattices are equal 
if we choose V = \JV s V d . Comparing now curves in both 
cases, we see that the data for the checkered lattice lie 
somewhat above the data for the striped lattice, joining 
together at densities 8 = 0, 0.5, and 1. Both models have 
the same V — > and V — > oo limits. 

In the Fig. 4 we plot the change in the shape of the 
density dependence of the diffusion coefficient as a func- 
tion of the ratio r. We keep V s = W s and V d = W d . Now 
the increase of the diffusion with decreasing parameter r 
does not lead to a singular behavior as r — > oo. Com- 
paring Fig. 3 with Fig. 4 we can see evident qualitative 
difference in the behavior of both systems. This differ- 
ence was not so clearly seen in Fig. 2, where results for 
different values of the rate V were plotted. 



V. MONTE CARLO SIMULATIONS 

In order to verify our analytical results for the be- 
havior of the diffusing gas on the nonhomogeneous two- 
dimensional surfaces we have simulated such systems 
by using Monte Carlo dynamical approach. We have 
used Boltzmann-Matano analysis of the shape of step- 
like density profile, after letting it decay via diffusion 
proces s 29 ' 30 . Results of this analysis are shown in Fig. 
5. Comparison for two different jump rates along y— 
direction are shown. These rates have been chosen in 



such a way that V = ^/V s V d , so the analytic results merge 
at 9 = 0.5. In the right panel of Fig. 5 we see system 
with the same jump rates in both x- and y-directions. 
The highest diffusion coefficient is two times larger than 
the lowest one. The same difference can be seen for the 
Monte Carlo data, even if noise of the results is large. If 
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FIG. 5: Comparison of Monte Carlo data with analytical 
curves for the diffusion coefficient D x dependence on the to- 
tal density of the system 9. Left panel shows the results for 
V s = 10 W s and right panel for V„ = W B - Simulation data for 
striped system are plotted by full circles and analytical curves 
by dashed line, whereas simulations for checkered lattices are 
plotted by open circles and analytical curves by full line. 



we increase rate of jumps along ^-direction by a factor of 
ten, then the expected ratio of the highest and the low- 
est diffusion coefficient increases to four(left panel). This 
trend is confirmed by simulations shown in the left panel. 
Unfortunately due to the high level of noise the difference 
in the behavior between the two types of analyzed lattices 
is not clearly visible in the simulation results. However, 
it is evident that the behavior of the diffusion coefficient 
as a function of the density changes with the increasing 
rate of jumps along y-axis in both cases in the manner 
consistent with our analytical results. 



VI. CONCLUSIONS 

We have shown that the recently formulated varia- 
tional approach to the collective diffusion is an effective 
and promising method of calculation of the diffusion co- 
efficient in two-dimensional systems. Here, we have used 
this method to describe behavior of a system of particles 
on a nonhomogeneous potential landscape. The resulting 
density dependence of the collective diffusion coefficient 
as a function of all rates that are present in the system 
is given by a simple analytic formula. This dependence 
agrees with Monte Carlo simulation results obtained for 
selected systems. We show that dynamical behavior of 
two-dimensional system is interesting and far from trivial 
even if site blocking is the only interaction that particles 
experience. It appears that in contrast to a single parti- 
cle system, collective diffusion in x- and y-directions of 
a system of particles depend on each other in the sense 
that the diffusion coefficient along one direction strongly 
depends on the rate jumps of particle in the direction 
perpendicular to it. This effect is induced by the site 
blocking, because it is not present in the system of many 
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independently moving particles We have demonstrated 
that the one-dimensional character of diffusion changes 
continuously when the rates of jumps are varied in di- 
rection perpendicular to the one along which diffusion is 
observed. We can understand this as a result of an acti- 
vation of alternative diffusion pathways, when the direct 
pathway is blocked. The transition from the one to two 
dimensional behavior is highly nontrivial even if particles 
do not interact. 
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